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^ ■ Abstract 
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In a recent paper Abramowicz and Kluzniak [1] have discussed the problem of 
epicyclic oscillations in Newton's and Einstein's dynamics and have shown that 
Newton's dynamics in a properly curved three-dimensional space is identical to test- 
body dynamics in the three-dimensional optical geometry of Schwarzschild space- 
time. One of the main results of this paper was the proof that different behaviour 
of radial epicyclic frequency and Keplerian frequency in Newtonian and General 
Relativistic regimes had purely geometric origin contrary to claims that nonlinearity 
of Einstein's theory was responsible for this effect. 

In this paper we obtain the same result from another perspective: by repre- 
senting these two distinct problems (Newtonian and Einstein's test body motion 
in central gravitational field) in a uniform way — as a geodesic motion. The so- 
lution of geodesic deviation equation reproduces the well known results concerning 
epicyclic frequencies and clearly demonstrates geometric origin of the difference be- 
tween Newtonian and Einstein's problems. 

Key Words: Geodesic deviation, Jacobi geometry, epicyclic oscillations, Schwarzschild 
spacetime 
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1 Introduction 



In a recent paper Abramowicz and Kluzniak [1] have discussed the problem of epicyclic 
oscillations in Newton's and Einstein's dynamics and have shown that Newton's dynamics 
in a properly curved three-dimensional space is identical to test-body dynamics in the 
three-dimensional optical geometry of Schwarzschild space-time. Their discussion was 
motivated by the theory of accretion disks around black holes and neutron stars which is 
based on assumption that accreting matter moves on nearly circular geodesic trajectories. 
One of the strong field effects that should be present in this context, as pointed out by 
Abramowicz and Kluzniak in a recent series of papers [2] , is the possibility of parametric 
resonance (preferably 3:2) between vertical and radial epicyclic frequencies of perturbed 
circular orbits. It has been conjectured [2] that this effect is indeed responsible for the 
observed double peaked QPOs [3]. 

It is quite well known that in the Schwarzschild (or Kerr) spacetime radial epicyclic 
frequency u r is lower than orbital frequency ujk (and vanishes at the marginally stable 
orbit) unlike in the Newtonian gravity where these two frequencies are equal. As recalled 
by the authors of [1] many people attributed this different behaviour to the nonlinearity 
of the Einstein's theory of gravity. Therefore one of principal motivations for [1] (besides 
making a brilliant use of the so called optical geometry [4]) was to demonstrate the purely 
geometric origin of this effect. In order to achieve this Abramowicz and Kluzniak have 
represented the Einstein equations (in optical geometry) for the motion on a circular 
orbit in Schwarzschild space-time in the form of Newton's equations in certain curved 3- 
dimensional space. Then they were able to calculate the epicyclic frequencies in a uniform 
way (i.e. from the same equation) and show explicitly that the aforementioned difference 
{uj r < ujk in Einstein's gravity vs. u r = uk in Newton's gravity) has purely geometric 
origin. 

In this paper I will obtain the same result from another perspective: by represent- 
ing these two distinct problems (Newtonian and Einstein's test body motion in central 
gravitational field) in a uniform way — as a geodesic motion. The difference in achieving 
"uniformity" is that whereas in [1] it was the same functional form of the equation in our 
case it will be the same geometric representation of the problem. 

2 Relativistic epicycles from the geodesic deviation 
equation 

Before going to the details let us start with some general comments. First of all, the 
problem of epicyclic frequencies has nothing to do with nonlinearity of Einstein's equations 
just because the Einstein's equations in general are dynamical equations for evolving the 
3-geometry (see e.g. [5]). In the problem of epicyclic oscillations around circular orbits 
one has a kinematic problem of test bodies moving in static spacetime — the geometry 
is static and defined a priori. Hence the relevant question is how do the adjacent orbits 
of test particles behave. 

The transition from Newtonian gravity to the Einstein's picture can be summarized 
in the following way. Newton's explanation why the planetary orbits are curved (circular, 
elliptical, parabolic or hyperbolic - for comets) was that it is the force of gravity from 
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the central body (the Sun) that makes them curved. In Newton's theory the nature of 
the force of gravity remained unexplained - it was taken for granted. Of course basic 
properties of the gravity force were explained e.g. the inverse square law, but not its 
nature. On the other hand, Einstein attempted at explaining the nature of gravity - there 
is no force field but the presence of massive central body makes the spacetime curved. 
The motion of test bodies takes place along geodesies; they are in a free motion but in a 
curved spacetime, that is why their trajectories are curved. 

Therefore the (general relativistic) problem of epicyclic orbital oscillations in Schwarzschild 
spacetime is exactly the problem of geodesic deviation in Schwarzschild geometry. Stable 
circular orbits are stable in the sense that geodesic deviation equation solved along such 
circular orbit has oscillatory solutions. 

One of the most recent papers presenting solution of the geodesic deviation equation 
for trajectories close to circular orbits in the Schwarzschild space-time is [6]. We will 
sketch main steps leading to the formula for radial epicylic frequency in Schwarzschild 
metric referring the interested reader to [6] for computational details. 

In a pseudoriemannian manifold with the line element 

ds 2 = g iiv dx ii dx v (1) 

the curve 7 S := x^(s) parametrized by the affine parameter s is a geodesic if the tangent 
vector = ^ is paralelly transported along 7 S : 

flr = -sr +r "»" v (2) 

Then consider close geodesic %. The vector ^ representing the separation between 
the geodesic 7 S and an adjacent geodesic % satisfies the geodesic deviation equation 

■^r = -W*r (3) 

Now, let us take the Schwarzschild metric 

dg2 = (1 _ 2GM )c 2 rft 2 _ 1 dr 2 - r 2 (d6 2 + sin 2 9d<p 2 ) (4) 

c r (1 c2j , ) 

It is well known that circular orbits r = R = const, are geodesies in the metric (4) and 
test particles move along such orbits with the angular velocity ujk (Keplerian frequency) 
given by the formula [5] 

-k = ^T (5) 

If one considers a nearby geodesic (with respect to the circular one) and asks how does the 
separation between these two behave, the answer would come from solving the geodesic 
deviation equation (3). Technically one should express the components of the tangent 
vector m m = (u f , u r , u e , u v ) as well as the components of the Riemann tensor R^ pa (also 
the Christoffel symbols while calculating a covariant derivative, etc.) as evaluated along 
the circular orbit. Then one gets the system of four second-order differential equations. 
One of them — for £ 9 component — reads = —oj 2 K ^ e and is decoupled from the rest. 
The remaining three form the system of coupled linear second order differential equations 
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with constant coefficients (detailed calculations can be found in [6]) and the characteristic 
equation for this system (written in matrix form) reads: 

leading to the value of radial epicyclic oscillations with the frequency: 

2 CM , 6GM 

In summary, the conclusion from solving the geodesic deviation equation in Schwarzschild 
spacetime is that behaviour of geodesies close to circular orbits can be represented as 
a superposition of (epicyclic) oscillations around circular orbit with two characteristic 
frequencies: the vertical epicyclic frequency - equal to Keplerian frequency ujk of the 
reference orbit and radial epicyclic frequency u r . 



3 Classical mechanics represented as problem of geodesies 

It is well known that variational principles of classical mechanics make it possible to 
formulate the dynamics of Hamiltonian systems as geodesic flows on some Riemannian 
manifold. This picture comes quite naturally from the Maupertuis-Jacobi least action 
principle (motion of the system with fixed energy E between q' and q" takes place along 
a path 7 minimizing the Maupertuis-Jacobi action): 

5S = 5 p ^E - V(q)Jaij dtfdqi = (8) 

Jq' 

and its formal resemblance to the variational formulation of geodesies in Riemannian 
geometry as curves extremalizing the distance. This is the simplest way to see desired 
correspondence and it is quoted in many textbooks on classical mechanics (e.g. in [7]). 
Below we give some steps along a straightforward "brute force" derivation of this result 
which could be instructive in seeing the role of time reparametrisation which is necessary 
in this picture. 

Consider the classical mechanical system described by the Hamiltonian 

H = H(p,q) = ^p tPj + V(q) (9) 

The equations of motion for the q variables with respect to time parameter t (in Newtonian 
physics one has an absolute time) and corresponding to the Hamilton equations may be 
written as 

q J +r\ s q s q k = -a^^ (10) 

where T J ks are the Christoffel symbols calculated with respect to metric and dots denote 
t - time derivative. Due to the force term this is, obviously, not a geodesic equation. It 
is simply the Newton's second law restated. The momentum variables are just linear 
combinations of velocities Pi = a^q 3 . Transformation to a geodesic motion (i.e. free 
motion in a curved space) is accomplished in two steps: (1) conformal transformation of 
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the metric a^-, and (2) change of the time parameter along the orbit. More explicitly we 
equip the configuration space with the metric - the so called Jacobi metric 

9ij = 2{E - V{q))a l3 (11) 

(note that is read off from the kinetic energy term in the Hamiltonian and (in general) 
is allowed to vary as a function of the configuration space variable, = (%((?)) )• Let 
us also call this Riemannian space i.e. configuration space accessible for the system and 
equipped with Jacobi metric — the Maupertuis- Jacobi manifold. 

With respect to the metric (11) and the time parameter t it is not easy to see that the 
orbits are geodesies since there is a term appearing on the right hand side of the equation, 



dT> q +T ^Jt q Jt q = -^W)Jt q W {q) dt 



where now denote the Christoffel symbols associated with the Jacobi metric. However, 
if we reparametrize the orbit q l = q l (s) in terms of the parameter s defined as 

ds = 2(E-V)dt (13) 

the orbits will become affinely parametrized geodesies, i.e. the configuration space vari- 
ables q l satisfy the well known geodesic equation 

^ ,+Tl 4/h k =° (14) 

with no force term on the right hand side. 

The information about the original force acting on the particle (as described by the 
potential V(q) in the Hamiltonian (9)) has been encoded entirely in the definition of the 
Jacobi metric (11) and the definition of the new parameter s in (13) parametrizing the 
orbit. 

Contemplating how nearby orbits behave (i.e. the local instability properties), it is 
natural to consider the geodesic deviation equation which describes the behavior of nearby 
geodesies (14). 

This can be derived in a usual manner by subtracting the equations for the geodesies 
q % (s) and ?*(s)+£ l (s) respectively or simply by disturbing the fiducial trajectory (pi(t), q l (t)) : 

Pi(t) = Pi(t) +Vi(t), 

q\t) = f{t)+C(t) (15) 

and substituting this directly into Hamilton's equations. In this way we also arrive mo- 
mentarily, though tediously, at the geodesic deviation equation for the separation vector 

£, 

^| = -iJ5„«W (16) 

Here u l = Dq l /Ds is the tangent vector to the geodesic, is the separation vector 
orthogonal to u. Note, that the covariant derivative D/Ds and Christoffel symbols are 
calculated with respect to the Jacobi metric (11). 
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4 Keplerian circular orbits represented as geodesies 
- epicyclic frequency in Newtonian's regime 

Geometric formulation of the Kepler problem is very simple. The two body problem in 
Newtonian gravity is essentially two-dimensional. Therefore the Jacobi metric (in polar 
coordinates) reads: 

ds 2 = 2{E-V{r)){dr 2 + r 2 d^ 2 ) (17) 

where: V(r) = — Moreover all information carried by the Riemann curvature tensor 
is captured by the Gaussian curvature. 

In this case the geodesic deviation equation (for an orthogonal Jacobi field £ of geodesic 
deviation g(£,u) = 0) reads: 

^- 2 C + K G e = (18) 

where Kq is Gaussian curvature of respective Maupertuis- Jacobi manifold. 

It is a simple exercise to calculate the Gaussian curvature of Maupertuis- Jacobi man- 
ifold for the Kepler problem. Let us denote f(r) 2 := 2(E — V(r)) so that Jacobi metric 
reads ds 2 = f(r) 2 (dr 2 + r 2 d(f 2 ), then let us consider the one-forms uj 1 := f(r) dr and uo 2 := 
f(r)r dip. Now, it is quite obvious that uj 1 is a closed form, and duo 2 = — j^^^-^-^lu 2 Alu 1 . 
Then from Cartan's equations one can easy read off the Gaussian curvature 

K= 1 d ( 1 d(rf(r)) \ _ EGM 

° f(r) 2 rdr\f(r) dr ) A(rE + GM) 3 { ' 

Note that since the kinetic energy is positive definite T := E — V(r) > the term 
rE + GM is also positive and the sign of Gaussian curvature is determined by the sign 
of the energy. It is negative for E > i.e. for hyperbolic orbits and positive for E < 
i.e. for bound motion. The meaning of this results is that, in the first case, scattering 
of test particles on the center has sensitive dependence on initial conditions (problem is 
equivalent to congruence of geodesies on negatively curved manifold) while in the second 
case (equivalent to circular or elliptic orbits) the disturbed trajectories execute Keplerian 
epicyclic oscillations around the orbit of reference. 

For the circular orbit of radius R (where E = 1/2V(R)) we have: 

K ° = am (20) 

and (by virtue of geodesic deviation equation) K G = oj 2 3 where cuo,j is the (Keplerian) 
epicyclic frequency in Jacobi geometry i.e. with respect to natural "time" s along the 
geodesic. 

Recalling that ds = 2(E — V(R))dt = ^jfdt one can easily recover Keplerian epicyclic 
frequency in Newtonian picture i.e. with respect to the Newtonian time t 



2 G M 2 
= = Uk (21) 

and it turns out to be equal to the orbital (Keplerian) frequency ujk- 

Let us also remark on the equation (20). The quantity in the denominator is in fact 
equal to the Keplerian angular momentum (per unit test body mass) squared. In our 
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approach it turns out that for circular orbits the specific angular momentum squared is 
equal to the inverse Gaussian curvature of the Jacobi manifold in which the motion takes 
place along geodesies. In the picture developed by Abramowicz and Kluzniak the specific 
angular momentum had an interpretation of the geometric mean of the gravitational 
radius of central body and the radius of curvature of particle's orbit. 

5 Conclusion 

Conclusion of this note is that by applying a uniform representation of the problem i.e. by 
representing physical trajectories as a problem of geodesies in some manifold one can see 
the geometric origin of the difference in epicyclic frequencies (describing the behaviour of 
trajectories adjacent to circular orbits) calculated in Newtonian and General Relativistic 
regimes. This result is in agreement with that of [1] although it has been derived in a 
different manner. One may say that instead of applying the Feynman's principle "the 
same equations have the same solutions" we have sucessfully applied the principle of 
"comparing comparable things" by working in the same geometric representation of the 
problem. 
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